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LASRY-LIONS, LAX-OLEINIK AND GENERALIZED 
CHARACTERISTICS 

CUI CHEN AND WEI CHENG 


Abstract. In the recent works [9] and m, an intrinsic approach of the prop¬ 
agation of singularities along the generalized characteristics was obtained, even 
in global case, by a procedure of sup-convolution with the kernel the funda¬ 
mental solutions of the associated Hamilton-Jacobi equations. 

In the present paper, we exploit the relations among Lasry-Lions regulariza¬ 
tion, Lax-Oleinik operators (or inf/sup-convolution) and generalized charac¬ 
teristics, which are discussed in the context of the variational setting of Tonelli 
Hamiltonian dynamics, such as Mather theory and weak KAM theory. 


1. Introduction 

Suppose H = H{x,p) : K” x K” ^ K is a Tonelli Hamiltonian, i.e., H is 
convex in p with superlinear growth condition. Let u : K" —K be a (global) 
viscosity solution of the Hamiltonia-Jacobi equation 

(1.1) H{x, Du{x)) = 0, a; G K". 

Such a solution u is locally semiconcave (with linear modulus) on R". We denote by 
D'^u{x) the superdifferential of m at a; (see, for instance, [13)1. which is a compact 
convex set in R”, and we call x G R" a singular point of u if D~^u{x) is not a 
singleton. Certain “singular dynamics” was interpreted in [2] by a Hamiltonian 
inclusion 

x(s) £ CO Hp{x{s), D'^u{x{s))), a.e.sG [0,t], 
and such a Lipschictz arc x is called a generalized characteristic. If xq is a singular 
point of u and 

(1.2) 0 0 CO Hp{xo, D^u{xo)), 

then the associated generalized characteristic x(t), t G [0, r], is composed of singular 
points of u. In the recent works [3] and [H], the propagation of singularities along 
generalized characteristics in m has been explained by an intrinsic variational 
approach (see, also, [2] [4] [TO] [13] [M] [32] for the approach from Control 

theory or PDE), which is motivated by Mather theory and weak KAM theory. 

Let us recall the aforementioned results in [S] and m at first. Let u G C(R"'), 
for any t > 0, Ti, the Lax-Oleinik operator of positive type, is defined as 

(1.3) ftu{x) := sup {u{y) - At{x, y)} := sup 'ijjtiy), x G R”, 
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where 


At{x,y)= inf / L{-f{s),'y{s))ds 

Jo 


with 

K,y = {7 € Vhi’i([0,t],R") : 7(0) = ai,7(t) = y}- 
Here L is an arbitrary Tonelli Lagrangian on R" with H its Fenchel-Legendre 
dual, and it is well known that m is also called sup-convolution or Lax-Oleinik 
opterators in the literature. In [9] and |12] . the authors shown that the maximizers 
in such a procedure of sup-convolution give exactly a local or global generalized 
characteristic starting from a singular points of u under suitable conditions. 

In the present paper, we will explain the connection between generalized char¬ 
acteristics and the well-known Lasry-Lions regularization at first. Throughout this 
paper, we suppose that L satisfies condition (LI) and (L2) (see Section 2). 

Let M be a closed manifold, t > 0 and u : M —)• R is any semiconcave 
function, the following properties are already known (see, for instance, [B], [22] 1 

(PI) Ttu belongs to class for 0 < t < to with to is a constant dependent on 
the constant of semiconcavity of u. 

(P2) TtU is decreasing on (0, -boo) if m is a viscosity subsolution of Hamilton-Jacobi 
equation 

H{x, Du{x)) = a{0), x£ M, 

where «(■) is Mather’s a-function. Moreover, TtU tends to u uniformly as 
t 0+. 


In this paper, we also have 


Theorem 1.1. Suppose u : R" —^ R is a semiconcave function. Then there exists 
0 < to ^ 1 such that 

(PS) Let xq € R" and L(xo,0) ^ 0, then Ttu(xo) is increasing on (0, -boo) and 
Ttu(xo) = u{xo). Consequently, if 

L{x, 0) <0, Vx S R", 

then TtU tends to u uniformly on any compact subset as t —> 0“''. 

(P4) Letx G R", suppose that the function ipt defined in (11.31) attains the maximizer 
yt in B(x, R(x,t)) for all t > 0 where R{x,t) > 0 is defined in (jA.bjl . then 
DTtu{x) = po, where po is the unique element with minimal energy: 

(ME) H{xo,p) ^ H{xo,Po), VpeII+M(xo). 

It is worth noting that the minimal energy condition in (IMEIl is the same as the 
inititial condition on the velocity of the generalized characteristic obtained by the 
intrinsic approach in [5] and m, see also Proposition 13.II 

In the rest part of this paper, we try to exploit the nature of the singularities 
of u by the procedure of inf-convolution. As pointed out in [9] , the inf-convolution 
defined by 

(1.4) Ttu{x) := ini {u{y) + At{y,x)}, a: e R”, 

yeR" 

is not the dual procedure of sup-convolution. But, it is still meaningful for study 
the critical points of the local barrier function 

(ptix) = u(x) + Atix,xo) 
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with respect to any fixed point xq. Recall that a point x G R" is a critical point 
of a locally semiconcave function rt if 0 G D'^u{x). Comparing to the local barrier 
function 

-iptix) = u{x) - At{xo,x), 

the function ■i/ij only admits a unique critical point (maximizer) for small time t > 0 
since the convexity properties of the fundamental solutions At{xo,x) (see, Appendix 

A). 

Along this line, given a singular point x of u, using a nonsmooth critical point 
theorem by Shi dsn), we obtain a critical point of the local semiconcave function 
cj)t, which is not a global minimizer of (j)t determined by classical characteristic 
passing to x. 

Theorem 1.2. Let u be a Lipsehitz viscosity solution of HH), t > 0, and let x G R" 
be a singular point of u. Suppose there exist finite many elements in D*u{x), say 
D*u{x) = {pi,... ,pk} with k ^ 2, then there exist critical points {xf} of (ft (not 
global minimizers) such that, for 1 ^ i,j ^ k, i ^ j, each critical point Xt = xf 
has the following dichotomy: 

(a) Xt is a differentiable point of (ft and there exists a local minimal curve connect¬ 
ing Xt and X. More precisely, there exists a curve 7 : (— 00 , t] —>■ R" such that 
7 ( 0 ) = Xt, 7 (t) = X and the restriction of ^ on (—oo,0] is a {u, L,0)-calibrated 
curve, hut 7 is not a {u, L,0)-calibrated curve on (—oo,t]; 

(b) Xt is a singular point of u. 

From the theorem above, the location of singularities affords possible information 
to construct “local” minimal orbits for Tonelli Lagrangian systems, which is totally 
unknown before. In the previous works of variational approach of Hamiltonian 
dynamical instability problem like Arnold diffusion (see, for instance, [7] , [19] , m, 
[El, US], m and [18jb the diffusion orbits shadow the variational minimizers 
which are not local ones. 

The paper is organized as follows: In section 2, we review some basic properties of 
viscosity solution of Hamilton-Jacobi equations. In section 3, we discuss the relation 
of the generalized characteristics give by the procedure of sup-convolution and 
Lasry-Lions regularization, then, we also discuss what happens for the procedure 
of inf-convolution. 

Acknowledgments This work was partially supported by the Natural Scientific 
Foundation of China (Grant No. 11271182 and No. 11471238), the National Basic 
Research Program of China (Grant No. 2013CB834100). The authors are grateful 
to Liang Jin for helpful discussions on the results of this paper. 

2. Viscosity solutions and weak KAM theory 

A function L : R" x R" ^ R is said to be a Tonelli Lagrangian if the following 
assumptions are satished. 

(LI) The Hessian ^:^{x,v) is positive definite for all {x,v) G R" x R". 

(L2) there exists a non-decreasing function 6 : [0,-|-oo) —>■ [0,-koo), 0{r)/r —> -koo 
as t ^ -koo, Co > 0 and ,Ci = Ci(x, i?) ^ 0 such that 

L{x, v) ^ 0(|u|) — Co, (x, v) G R" X R”, 

and 


\Lx{y,v)\ -k \Ly{y,v)\ ^ ci{x,R)e{\v\), {y,v) G B{x,R) x R”. 
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Let H : R” x K" —>■ K be the associated Tonelli Hamiltonian, i.e., 
H{x,p)= sup{{p,v) - L{x,v)}. 


Throughout this paper, we suppose L is a Tonelli Lagrangian with conditions 
(LI) and (L2). 

2.1. semiconcave functions. Let C R." be a convex open set, a function u : 

—>■ R is semiconcave (with linear modulus) if there exists a constant C > 0 such 
that 


(2.1) Xu{x) + (1 - X)u(y) - u{Xx + (1 - X)y) < ^A(l - A)|a; - yf 


for any x,y € and A G [0,1]. Any constant C that satisfies the above inequality 
is called a semiconcavity constant for u in 12. A function it : 12 —>■ R is said to be 
semiconvex (with linear modulus) if —it is semiconcave. A function it : 12 —> R is 
said to be locally semiconcave (resp. locally semiconvex) if for each x € 12, there 
exists an open ball B(x,r) C 12 such that it is a semiconcave (resp. semiconvex) 
function on B(x,r). 

Definition 2.1. Let it : 12 C R" —>■ R be a continuous function. We recall that, for 
any x € 12, the closed convex sets 



are called the (Dini) subdifferential and superdifferential of it at x, respectively. 

Definition 2.2. Let it : 12 R be locally Lipschitz. We recall that a vector p G R" 
is called a limiting differential of it at x if there exists a sequence {xn} C 12 \ {x} 
such that It is differentiable at Xk for each fc S N, and 


lim Xfc = X and lim Du{xk) = P- 



The set of all limiting differentials of it at x is denoted by D*u{x). 

The fundamental properties of the superdifferential of a semiconcave function 
are listed in the following proposition. The monograph [13] is a good reference for 
the topic of semiconcave functions and beyond. 

Proposition 2.3. Let u : LI G R" M. be a semiconcave function and let x G 12. 
Then the following properties hold. 

(a) D^u{x) is a nonempty compact convex set in R" and D*u{x) C 9Z2+it(x), 
where dD^u{x) denotes the topological boundary of D^u{x). 

(b) The set-valued function x D'^u(x) is upper semicontinuous. 

(c) If D'^u{x) is a singleton, then u is differentiable at x. Moreover, if D'^u(x) is 
a singleton for every point in 12, then u G C^(12). 

(d) Z2+it(x) = coZ2*it(x). 

(e) D*u{x) = { limi^ooPi : Pi G D^u{xi), Xi —^ x, diam (Z2+it(xi)) O}. 

From proximal analysis point of view, the following result characterizes the semi¬ 
concavity of a continuous function and its superdifferential. 
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Proposition 2.4. Let m : £7 —>■ R fee a continuous function. If there exists a 
constant C > 0 such that, for any x € LI, there exists p G R" such that 

C 

(2.2) u{y)ifu{x) + {p,y-x) + —\y-x\^, ^y € LI, 

then u is semiconcave with constant C and p G D'^u{x). Conversely, if u is semi¬ 
concave in LI with constant C, then (1221) holds for any x G LI and p G D^u{x). 

Finally, we introduce the concept of singularity of a semiconcave function. A 
point X G n is called a singular point of u if D^u{x) is not a singleton. The set of 
all singular points of u, also called the singular set of u, is denoted by E„. 

2.2. Fundamental solutions and viscosity solutions. Given x,y G R", we 

define 

= {7 e W^’\[0, i], R”) : 7(0) = x, j{t) = y} 

Let t > 0, we denote 

(2.3) Atix,y)= inf f L{'y{s),^{s))ds. 

Jo 

It is well know that the infimum can be achieved by (7^ curves. In the literature of 
PDFs, At{x, y) is called a fundamental solution of (12.41) . see, for instance, [29] . 

Throughout this section, we suppose the Tonelli Lagrangian L satisfies con¬ 
dition (L1)-(L2). We discussed the associated Nagumo type conditions and the 
essential regularity results of the fundamental solutions in Appendix |A| For the 
main regularity results we will use, see, Proposition IA.3I and Proposition IA.4I 
Suppose iJ is a Tonelli Hamiltonian, throughout this paper we will be concerned 
with the Hamilton-Jacobi equation 

(2.4) H{x,Du{x)) =0, xGR”. 

We recall that a continuous function u is called a viscosity subsolution of equation 
(1^ if, for any x G R", 

(2.5) H{x,p)^0, \/pGD^u{x). 

Similarly, u is a viscosity supersolution of equation (12.41) if, for any x G R” , 

(2.6) H{x,p)^0, ypGD~u{x). 

Finally, u is called a viscosity solution of equation (EH), if it is both a viscosity 
subsolution and a supersolution. 

Proposition 2.5. Any viscosity solution of the Hamilton-Jacobi equation (12.41) is 
locally semiconcave with linear modulus. 

Proposition 2.6. Exti7+u(x) = D*u{xj^for any viscosity solution u of (12.41) and 
any x G R". 

Proposition 2.7. Let x G R" and u : R” —>■ R fee a viscosity solution of the 
Hamilton-Jacobi equation (EH). Then p G D*u(x) if and only if there exists a 
unique curve 7 : (—oo,0] ^ R" with 7(0) = x which is which is a (u,L,0)- 
calibrated curv^ and p = Ly{x,'y{0)). 


^For any convex closed subset of R", we denote by ExtC the set of extremal points of C. 
^For the concept of dominated functions and calibrated curves, see, for instance, m 
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2.3. Generalized characteristic. The construction of the singular set or cut loci 
of viscosity solutions is a very important and hard problem in many fields such as 
Riemannian geometry, optimal control, classical mechanics, etc.. It is known that 
the study of propagation of singularities can go back to [1] for general semiconcave 
functions by the method from nonsmooth analysis. Some dynamical nature of the 
singularity was found by the concept of generalized characteristic. 

Definition 2.8. A Lipschitz arc x : [0,t] ^ R." is said to be a generalized charac¬ 
teristic of the Hamilton-Jacobi equation (12.4|) if x satisfies the differential inclusion 

(2.7) x(s) G coi7p(x(s),I?'''u(x(s))), a.e. s G [0,r]. 

A basic criterion for the propagation of singularities along generalized charac¬ 
teristic was given in [5] (see |141132) for an improved version and simplified proof). 

Proposition 2.9 ([2]). Let u be a viscosity solution of Hamilton-Jacobi equation 
(1^ and let xg G M". Then there exists a generalized characteristic x : [0,t] —)■ R." 
with initial point x(0) = xg. Moreover, if xg G Su, then x(s) G Su for all s G [0,t]. 
Furthermore, if 

(2.8) Q ^ coHp{xg,D^u{xg)), 
then x(-) is injective for every s G [0,t]. 

3. Procedure of sup-convolution and generalized characteristics 

Let be a Tonelli Hamiltonian on R". Recall the Lax-Oleinik operators Tt and 
Tt, i.e., for any u G C(R”), 

(3.1) ftu{x):= sup {u{y) - At{x,y))}, 

(3.2) Ttu{x) := ini {u{y) + At{y,x))}. 

y^R^ 

When taking H{p) = |pp/2 and the kernel At{x, y) = A\x — yp, the two operators 
above are closely linked to the so-called Lasry-Lions regularization procedure ([IS]) 
which is written in the form of sup-convolution and inf-convolution, respectively. 
This type of regularization is also called Moreau-Yosida regularization in convex 
analysis. A more detailed formulation can be found in [5] with respect to the 
aforementioned quadratic kernel. 

3.1. Procedure of sup-convolution and generalized characteristics. Re¬ 
cently, in |5] and the authors studied the intrinsic relation of propagation 
of singularities along the generalized characteristics and the following procedure of 
sup-convolution. 

Fix X G R", 0 < t ^ tg 1, then there exists R{x, t) > 0 such that, the function 

(3.3) iftiy) ■■= u{y) - At{x,y), y G B{x, R{x,t)), 

has a unique maximizer for each t G (0, to]) where At{x, y) is a fundamental solution 
with respect to the associated Tonelli Lagrangian L. 

Suppose that u{-) is semiconcave while At{x, •) is locally semiconcave (Proposi¬ 
tion ??) and convex when t G (0,to] iProposition IA.3II . say Ci > 0 (resp. —C'2(t)) 
is the semiconcavity (resp. semiconvexity) constant of u(-) (resp. At(x,-)). Note 
that, by Proposition IA.41 the constant C 2 (t) = '!/'t(-) is strictly concave 
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in B(x, R(x,t)) and consequently we have a unique maximizer for each t G (0,fo], 
which is also a unique critical point of tpt- 
Let us define the arc y : [0, to] —t K" by 


(3.4) 


y(t) 


X, t = 0, 
yt, tG(0,to]- 


If : [0,t] -G K" is the unique minimizer in the definition of At{x,y), we dehne 
(3.5) pt{s) := s G [0,t], 


the associated dual arc with respect to ^t(s). 


Proposition 3.1 ([9]). Let u be a locally semiconcave function and x G Su, the 
singular set of u. There exists ti ^ to such that the arc y : [0,ti] -G K" defined 
in (1^ is a generalized characteristic composed of singular points of u, i.e., y : 
[0,ti] —>■ R” is Lipschitz continuous, y(t) G for all t G [0,ti], and satisfies 

(3.6) y(T) G coiJp(y(T),D+u(y(r))), a.e. r G [0,ti]. 

Moreover, 

(3.7) y+{Q)=Hp{x,po), 
where po is the unique element of minimal energy: 

H{x,p) ^ H{x,po), VpGD'^u{x). 


3.2. Lasry-Lions regularization. In this section, we will explain the connection 
between Lasry-Lions regularization ('|25)1 and generalized characteristics first found 
in [5]. We only concentrate to the case of sup-convolution Ttu with u a locally 
semiconcave function. 

For t > 0, recalling that 

(3.8) ftu{x) := sup {u{y) - At{x,y))}, 

yGR" 

where u : K” —>■ R is any locally semiconcave function, and At{x,y) is the funda¬ 
mental solution with respect to any Tonelli Lagrangian L. 


Theorem 3.2. Suppose u : R” ^ R 6e a semiconcave function with constant C. 
Then there exists 0 < to 1 such that ifTtU is defined as in p.8l) . we have 

(PS) Fix X G R”, then Ttu{x) is increasing on (0, -|-oo) and limj_,,o+ Ttu(x) = u(x) 
ifL(x,0) ^ 0. Consequently, if 

L(x,0)<0, Va:GR’", 

then TtU tends to u uniformly on any compact subset as t —>• O'*'. 

(P4) Letx G R", suppose that the function ipt defined in (13.3|) attains the maximizer 
yt in B(x, R{x, t)) for all t > 0, then limt_,.o+ DTtu{x) = pq, where po is the 
unique element with minimal energy, i.e., 

(3.9) H{x,po) = min H{x,p). 

p^D'^u{x) 

(P5) In particular, when L has the form 

L{x,v) = ^{Av,v), xGR",uGR", 
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where A is an n x n symmetric and positive definite matrix. If t ^ kC~^, 
then, the functions u and Ttu have the same critical points and critical values 
where k > 0 is the smallest eigenvalue of A. 

Remark 3.3. The properties (PI) and (P2) (see the introduction) is already known 
(in the case of compact manifolds), see, for instance, [6] or Fathi’s book [22]. Since 
this is a local result, it is not hard to generalize to the manifolds using local charts. 
We collect the known results here just for the comparison interests like (P2) (in the 
introduction) and (P3). The property (P5) is a slight generalization of a known 
result (|5]). 

Remark 3.4. It is worth noting that the assumption in (P4) that the function tpt 
defined in (13.31) attains the maximizer yt in B{x, R(x, t)) for all t > 0, is not easy to 
be checked in general. Fortunately, if we consider a certain type of nearly integrable 
systems or mechanical systems, this condition holds. The readers can refer to [12j . 


Proof. Let x,y € K" and t > 0, for any 0 < s < t, it is easily checked that 
At{x, y) ^ As{x, y) + At-s{x, x). 


Taking the constant curve 7 (t) = x, t € [0, t — s], we have 

At-s{x,x)i^ / L(7(t),7(t)) dr = (t - s)L(a;,0). 

Jo 

Therefore, for any fixed x € R”, we have At{x, ■) < As{x, •) since L{x, 0) ^ 0, and 
thus, tpt{y) ^ i’siy) for all y G R”. This leads to the conclusion that Tsu{x) ^ 
Ttu(x) if 0 < s < t. The uniform convergence result is a direct consequence of 
Dini’s Lemma on monotone sequence of continuous functions. This completes the 
proof of (P3). 

Now, we turn to the proof of (P4). Fix x G R" and t G (0, to]- Adopting the same 
terminology as before, since 'i/’t(') attains the maximum at j/ = yt G B{x, R{x,t)) 
and ^t G F* is the minimal curve in the definition of At{x, yt), we have 

Lv{^t{t),it{t)) = DyAt[x,yt) G D+u{yt), 

since the results in Proposition |A3| and 0 G D^iptiyt)- Moreover, we have that the 
family {^t}te(o.to] is equi-Lipschitz, by ProDOsition lA.2l 
Therefore, we have 


itit) -X 


-6(0) 


< 7 / l6(s) -6(0)1 ds 


t 

c 


^ — / s ds = — t. 

2 


Thus, we obtain 

Vo = lim vt = lim ^t(O) = lim it{t), 

t->0+ t->.0+ t->0+ 

where Vt = {yt — x)/t. Since m is a locally semiconcave function, thus, by the 
monotone property of semiconcave functions (see, e.g., [13]), we have 


(3.10) {p-Ly{^t{t),it{t)),vt)+ tC\vt\'^\/pGD+u{x). 


Taking limit in (I3.10L then 

(3.11) {p,vo) ^ {Ly{x,vo),vo), Vp G D'^u{x). 
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In other words, 

(3.12) H{x,p) ^ {Ly{x,Vo),VQ) - L{x,Vo) = H{x,po), \lp&D^u{x), 

where po = Ly(x, vq) € D^u{x), by the upper semicontinuity of the set valued func¬ 
tion X D^u{x), is the unique element solve the associated optimization problem 
(I3.12|) . and limt_,.o+ DTtu{x) = limt_>o+ L„(^t(0),^t(O)) = po- This completes the 
proof of (P4). 

For the proof of (P5), note that, in our case, the minimal curve ^t(s) = ■ s, 

thus, by (13.101) . we have 

{p - Avt.vt) + tC\vt\^ ^ 0, \/p€D^u{x). 

If 0 G D'^u{x), take p = 0 in the inequality above, then it follows 

{-KVt,vt) + fClutp ^ {-Avt,vt) + tc\vt\^ ^0, Vp G D+u{x), 

where k > 0 is the smallest eigenvalue of A. Therefore, 

(tC — K)|Ttp > 0, 'ip&D'^u{x). 

If t ^ kC~^, then vt = 0, yt = x and u{x) = Ttu{x). Conversely, if 0 = DTtu{x), 
then Vt = 0 and yt = x. It follows 0 G D^u{x) which proves (P5). □ 

3.3. What happens for the inf-convolution. In this section, we will discuss 
the procedure of inf-convolution, that is, let rt be a locally semiconcave function on 
M", and let L be a Tonelli Lagrangian, for any fixed x G K.", define 

(t>t{y)-=u{y)AAt{y,x), p G K". 

It is worth noting that </>( is the sum of two locally semiconcave functions, and it is 
also locally semiconcave consequently. 

For the convenience of our discussion, we suppose that m is a global viscosity 
solution of the Hamilton-Jacobi equation 

(3.13) H{x, Du{x)) = 0, xG R", 

where H is the associated Hamiltonian with respect to L. 

At this stage, we have 

u{x) = Ttu{x) = inf (j>t{y) 

yeR" 

for alH > 0 by well known facts from weak KAM theory. 

Lemma 3.5. Let u be a viscosity solution of 1)3.131) . and (pt is defined as above for 
t > 0. Then for t > 0, there exists Zt such that 

Mzt)= ini My)- 

yeR" 

Proof. This is actually obvious. Indeed, by Proposition 12.71 for any t > 0, and 
p G D*u{x), there exists a curve 7 : (—oo,t] —>■ R" such that 'y{t) = x, p = 

Lv(nit),fi{t)) and 

u(7(t)) - u(7(s)) = [ L(7(T),7(r)) dr, Vs < t. 


Take zt = 7(0), then we have the expected result. 


□ 
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Now, we can impose such a question: Is the aforementioned procedure of inf- 
convolution efficient for tracking the information of the propagation of singularities 
along generalized characteristics? 

We will try to answer this question using the technique from nonsmooth critical 
point theory , see also m for the applications by standard using Lasry-Lions 
regularization. 

Lemma 3.6. Let u be a viscosity solution of (|3.13|) and the function (ft be defined as 
above for any fixed x € K” and t > 0. Then there exists a one-to-one correspondence 
between p € D*u(x) and the global minimizers Zt of (ft for all t > 0. 


Proof. Let Zt € K” be a minimizer oi (ft, t > 0, then (ft is differentiable at Zt since 
(ft is locally semiconcave. Thus, zt is a differentiable point for both u and At{-,x). 
Consequently, there exists two curves 71 : (—oo,0] —K” and 72 G such 
that 


7i(0) = 72(0) = Zt , 
p = Du{zt) = Li,(7i(0),7i(0)), 
p' = DxAt{zt,x) = -L„(72(0),72(0)), 

by Proposition 12.71 and Proposition IA.41 and p -\- p' = 0 since Zt is a critical point 
of (ft- Moreover, 71 is a {u, L, 0)-calibrated curve, i.e., for any s > 0, 

w(7i(0)) -w(7i(-s)) = f L(7i(r),7i(r)) dr, 


and, similarly, 

u{x) - 11 ( 72 ( 0 )) = At( 72 ( 0 ),a:) = f L( 72 (r), 72 (t)) dr. 

Jo 

By the juxtaposition of 71 and 72 , we define 

r,t(T) = l 

^ I 72 (r), 0 < T < t. 

It is clear that rjt is a curve on (—oo,f] with ptlt) = x, and 

u{x) - u{r]t{-s)) = / Lfqt{T),f]t{T)) dr, s >0 
J — S 


which follows that r]t is also a (m, L, 0)-calibrated curve, and such a {u,L,0)- 
calibrated curve passing through zt with x the terminal datum is unique. Therefore, 
the correspondence between zt and pt is one-to-one. 

The rest of the proof is a direct consequence of Proposition [5771 D 


Now, we fix a point x € K". 

(1) If X is a differentiable (or regular) point of u, then D*u{x) = {Du{x)}, and (ft 
has a unique global minimizer zt which determines a unique (u, L, 0)-calibrated 
curve passing though zt with x the terminal endpoint point. 

(2) If X is singular point of u, it become relatively complicated. Let 

Zx,E = {p G M" : H{x,p) < E}, 

which is a non-empty compact and convex set when the energy E, say E = 0, 
is suitably chosen. It is known that D*u{x) = TAit D'^u{x), the set of extremal 
points of D'^u{x), by Proposition 12.61 This means the elements of D*u{x) is 
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exactly the set ExtD^u{x) which is located in the energy hypersurface dZx,E 
since is strictly convex. 

In the spirit of Lemnia l3.61 we want to look for the critical points of </>(. A point 
X € R” is a critical point of a locally semiconcave function m if 0 £ D^u{x). To 
find the critical points of (pt besides the global minimizers as in Lemma 13.61 we can 
not apply the standard Lasry-Lions regularization directly since such a function 
(f>t is only locally semiconcave. Fortunately, recall a well known nonsmooth critical 
point theorem, see, for instance, m- We only need the result in the following finite 
dimension setting. 


Proposition 3.7. Let f : R." —>• R 6 e a locally Lipschitz function. Suppose that 
xi,X 2 £ R", X 2 ^ B(xi,r) with r > 0 such that 


max{f{xi),f{x 2 )}<bo< inf /, 

dB{xi,r) 


and define 


b = inf max fH{t)), 
7ertG[o.i] 


where F = {7 £ (^([0,1],R") : 7(0) = xi,7(l) = X 2 }. If f is coercive, then 
there exists X 3 such that /(ts) = b and 0 £ df^x^), where df(xz) is the Clarke’s 
generalized gradient of f at X 3 . 


The readers can refer to [2T] for the definition and properties of Clarke’s gener¬ 
alized gradients. Applying Proposition 13.71 to f = (ft above, we obtain 


Lemma 3.8. Let u be a Lipschitz viscosity solution of (13.131) . < > 0, and let x £ R" 
be a singular point of u. Suppose there exist finite many elements in D*u{x), say 
D*u(x) = {pi,... ,pk\ with k ^ 2, then there exist k distinct global minimizers 

Moreover, if 

(3.14) bij = inf max (^(( 7 ^ ( 5 )), 1 i, j if k, i ^ j, 

-yijeTij sG[o,i] 

where Ttj = {7 £ (^([0,1],R") : 7 ( 0 ) = = z{}, then, for each pair of {i,j) 

with i ^ j, there exists a third critical point x’f of (ft such that (ft{pcf) = bij > 
minygRr. (ftiy). 


Proof. We suppose k = 2 and the proof in the general case is definitely similar. 
Suppose {pi,P 2 } = D*u{x), and t > 0, then there exists two (m, L, 0)-calibrated 
curves rjl and rj^, and two global minimizers zf and z^ of (ft such that 

vl{t) = vUt) = X, ( 0 ) = z}, 772 ( 0 ) = z"^, 

by Lemma 13.61 Since z }, zf are two isolated (global) minimizers of local Lipschitz 
function (ft which is coercive by the superlinear growth condition on L, then a third 
critical point is obtained in the context of mountain pass method as in Proposition 
13.71 The rest part of the proof is a direct consequence of Proposition 13.71 and 
the facts the Clarke’s generalized gradient d(ft{-) coincides with the (proximal) 
superdifferential D^(ft{-) (see [l3j) since (ft is locally semiconcave. □ 


Therefore, we claim that if a; £ R" is a singular point of viscosity solution u, 
there exists a third critical point xt = x\‘^ of (ft determined by two isolated global 
minimizers z} and zf and Lemma 13.81 
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Theorem 3.9. Let u he a Lipschitz viscosity solution of (13.131) . < > 0, and let 
X G R" he a singular point of u. Suppose there exist finite many elements in 
D*u(x), say D*u{x) = {pi,... ,pk\ with k ^ 2, then there exist critical points 
{xl^} of (ft (not global minimizers) such that, for 1 ^ i, j ^ k, i ^ j, each critical 
point Xt = x^^ has the following dichotomy: 

(a) Xt is a differentiable point of ft nnd there exists a local minimal curve connect¬ 
ing Xt and x. More precisely, there exists a curve 7 : (—00, t] —>■ R" such that 
7(0) = Xt, 7(i) = X and the restriction of ^ on (—oo, 0 ] is a (u, L,0)-calibrated 
curve, but 7 is not a {u, L,0)-calibrated curve on (—oo,t]; 

(b) Xt is a singular point of u. 

Proof. The existence of such critical points of ft is a direct consequence 

of Lemma l3.81 The critical points {x]^} are not global minimizers of ft since each 
global minimizer zl, 1 ^ i ^ k, is isolated. 

Suppose Xt is not a singular point of ft, thus both u(-) and At{-,x) is differen¬ 
tiable at Xt and 

Du{xt) D„;At{xt, x) = 0 

since Xt is a critical point of ft- Let p = Du{xt) and p' = —DxAt{xt, x), then there 
exist two curves 71 : (—00,0] —)■ R" and 72 : [0, t] —)■ M" such that 

7i(0) = 72(0) = a;*, 72(0 = a:, 

and 

P = L«(7i(0),7i(0)) = - p ' = Lt,(72(0),72(0)). 

It follows that 7i(0) = 72(0) and 7, the juxtaposition of 71 and 72, is a curve 
which is an extremal. But, 7 : (—oo,t] —^ R is not a (u, L, 0)-calibrated curve, 
otherwise, 7(0) = x* is a global minimizer of ft by Lemma 13.61 □ 

Remark 3.10. It is not clear, even when f > 0 sufficient small, if the critical point 
x'f of ft found by Lemma 13.81 is close to the singular point x of u. It is not hard 
to prove that when the positive time t tends to 0, the global minimizers zl and zf 
tend to X along the direction determined by the associated limiting differentials in 
D*u{x) respectively. We hope to dig out more information from this approach in 
the future. 

Appendix A. Regularity properties of fundamental solutions 

For the details of the proofs of the results in this appendix, the readers can refer 
to [9] and [12], or |8] under certain special conditions. 

Proposition A.l. Let 0<t^l, i?>0 and suppose L satisfies condition (LI) 
and (L2). Let ^ € F* be a minimizer for At{x,y), x € R", y G B(x,R), and let 
p{s) be the dual arc off{s). Then we have 

sup \i{s)\ if A{x,R/t), sup |p(s)| ^ A(x, sup \^{s)\ sf A{x, R/t), 

sGfO.t] se[07] sGfO.t] 

where A(x, •) is non-decreasing and continuous. 

Proof. For any t > 0, i? > 0, let x € R”, y € B{x, R) and ^ G ^l:,y ^ minimizer 

for At{x,y), i.e., At{x,y) = Jq L(((s), ((s)) ds. Denoting by ct € Fj, ,^ the straight 
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line segment defined by a(s) = x + |(y — x), s G [0,f], then by the Nagumo 
conditions in (L2) we have 


9{\i{s)\) ds - Cot ^ / L{C{s),iis)) ds ^ / L{a{s),&{s)) ds 


/o 


'0 


/o 


= / L[x+^-{y-x),^ 


JO 

j , s s y-x 




ds 




^cit|a;|0 


y-x 


t 

y-x 


t 


+1 max 

yGB(x,R),sG[ 0 ,t] 

+ tM{t,R) 


r / s .y-x 


^cit\x\9(R/t) + t max \L(x,v)\ := Ci{t, R). 

\x\,\v\^R/t 

By condition (L2), we have 

\L{x,v) - L{0,0)\ ^\L{x,v) - L{x,0)\ + \L{x,0) - L{0,0)\ 

^ci9{R/t)\v\ + C2(x). 

Thus, 

(A.l) Ci{t,R) ^ co{x,R)tKi{R/t) 

with Ki(s) = 0(s)(l + s) + 1. By the superlinear growth condition of 9, we have 
that ^ 

/ l^('S)| ds ^ C4{x,R)tK2iR/t), 

Jo 

where K 2 (s) = 1 + ki(s). Hence 

(A.2) |^(s)— [ \i(s)\ ds ^ C 4 {x, R)tK 2 iR/t), se[0,t], 


and 

(A.3) 


inf |^(s)| ^ / |C(s)| ds ^ C4{x,R)tK2{R/t). 

se[o.t] t Jo 


Now, we turn to estimate supj,gjQ |C(s)|. By condition (L2), we have 

^(le(s)l) ^L{^{s)A{s)) < i(e(s),0) + {L,{^{s)A{s))Ai.s)) 

=msl 0) + ( r L.iar), iir)) dr + L.(C(0), ^(0)), e(s) 


Note that we use the Euler-Lagrange equation in the last equality. By condition 
(L2) and the estimates above, we have 

[ \Lxi^iT),i{T))\ dr [ ci9{\i{T)\) dr ^ C4(x,R)tK2iR/t). 

Jo Jo 

For any 5 € [0,i], we also have 

l^«(C(0),C(0))| < \Ly{^{s),i{s))\ + [ |Ta:(C(T),C(T))| dr. 

Jo 

Then, by condition (L2) and (IA.3I) . it follows 

l-^«(^(0),C(0))| < ci9{\i{s)\) + C4{t,R)tK2{R/t), 
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and this implies 

|i«(?(0),C(0))l ^ci6»( ini \^{s)\)+C4{t,R)tK2iR/t) 

se[o,t] 

^Ci9{c3{x, R)tKi{R/t) + Cot) + C 4 {x, R)tK 2 {R/t). 

It follows there exists M > 0 and /r > 0 such that 

\i{s)\ < ^{ci6>(c3(a;,i?)tKi(i?/t) + Cot) + 2c4(t,i?)tK2(i?/t) +/i}. 

So, if t ^ 1, we have 

sup |^(s)| ^C' 2 (a;,it/t) 
se[o,t] 

As for the dual arc p(-), by (L2), we have 

sup b(s)| = sup |L„(^(s),^(s))| < C' 3 (a;,it/t). 

sG[0,t] se[o,t] 

We complete the proof by defining A{x, R/t) = max{(72{a;, i?/t), C' 3 (x, R/t)}. □ 

Fix a; e M" and suppose i? > 0 and L satisfies condition (L1)-(L2). In this 
case, the following observation is one of the key points of the results on the local 
regularity properties of At{x, y).. For any t > 0 and y e B{x,R), let ^t,y S F^ ^ be 
a minimizer for At{x,y), and pt^y be its dual arc, then we have 

sup |Ct.y(s)K A(a;,i?/t), sup |pt,y(s)K A(a;, i?/t), 
se[o,t] se[o,t] 

by Proposition lA.il Now, define 

K,, := B{x, A{x, 1)) X .8(0, A(a;, 1)) C K" x M”, 

K; := B{x, A{x, 1)) X .8(0, A(a;, 1)) C K" x (M")A 

Then, by defining a function R{x, •) : K" x (0,1] (0, oo), R{x,t) = we have 

(A.6) A(a;, 1/2) ^ A(x, 1). 

because of the monotonicity properties of A(x, •) and the continuity. So, if y G 
B{x, R{x,t)), and G i'x,y is a minimizer in the definition of At{x,y), then 

{^(s),p(s))}se[o.t].te(o.i] C K*, {^(s),^(s))}sg[o,t],te(o,i] C K*, 

Proposition A.2 ([H]). Fix any x G K." and t > 0 with R{x,t) defined as in 
(IA.6I) . If yt is the unique maximizer of-ipt in B{x,R{x,t)) for all t G (0,to]: o-nd 
f,t G F^ j^^ is a minimal curve in the definition of At{x,yt), t G (0,to]) then the 
family is equi-Lipschitz. 

The proof of the following result is similar to those in [T^] since the estimates 
involving certain first and second order partial derivatives of L or i? which are 
bounded on the a priori compact sets K* or K^,. The difference between the cases 
here and what in m is that the bound for the minimal curves and the dual arc is 
independent of x in the latter. 

Proposition A.3. Suppose L is a Tonelli Lagrangian satisfying (L1)-(L2). Fix 
any x G M", then there exists to > 0, such that for 0 < t ^ to, {t,y) M- At{x,y) is 
locally convex in 

S{x,to) = {(t,y) G M X R" : 0 < t ^ to, |y — a^l ^ R{x,t)}, 
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with R{x,t) defined in (IA.6I1 . 

More precisely, there exists constants Ci,C 2 > 0 such that, if y G B{x, R{x,t)), 
then, for \h\ <C 1 and \z\ <C 1, we have 

(A.7) At+h(x,y + z) + At-h(x,y-z)-2At(x,y) + 

Proposition A.4. Suppose L is a Tonelli Lagrangian satisfying (L1)-(L2). For 
any x € R", there exists tg > 0, such that the functions w : {t,y) i—>• At{x,y) and 
{t, y) I—!> Afiy, x) are both of class in 

S{x,to) = {{t,y) G M. X : 0 < t to, \y — x\ ^ R{x,t)}, 
with R{x,t) defined in (jA.6|) . for 0 < t ^ to. In Particular, for any t € (0,to]; 

(A.8) DyAt{x,y) =Ly{^{t),i{t)), 

(A.9) D,,At{x,y)=-L^{m,m). 

(A.10) DtAt{x, y) =- 

where ^ G is the unique minimizer for At{x,y) and Et^x,y is the energy of the 
Hamiltonian trajectory {^{s),p{s)) with p{s) = Ly{^{s),^{s)). 
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